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One of the most prominent features of the 7 -ray sky is the emission from our own Galaxy. 
The Galactic plane has been observed by Fermi-hAI in GeV and H.E.S.S. in TeV light. Fermi 
has modeled the Galactic emission as the sum of a complex ‘diffuse’ emission model with the 
predominately point source catalogs of IFHL and 2FGL, while H.E.S.S. has primarily detected 
extended TeV sources. At GeV energies, Galactic diffuse emission dominates the 7 -ray Milky 
Way but, as sources have hard spectra, it is likely their emission dominates at TeV energies. 
Generally the spatial shape and fraction of source emission compared to diffuse emission in the 
Galactic plane is not well known and is dependent on the source detection method, threshold and 
diffuse emission modeling methods used. 

We present a simple image-analysis based method applied to Fermi-LAI data from 10 GeV to 
500 GeV, covering a region of -e/- 5 degrees in Galactic latitude and -e/- 100 degrees in Galactic 
longitude, to separate source and diffuse emission. This method involves elongated filter smooth¬ 
ing, combined with significance clipping to exclude sources. We test the method against models 
based on the IFHL catalog and very simple model Galaxies to evaluate the response for an input 
of known fraction and shape of diffuse and source emission. 
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1. Introduction 

In this work, we study the separation of sources from Galactic diffuse emission using image- 
based techniques. These are applied to data at energies between 10 and 500 GeV in order to 
understand appropriate methods for automatic diffuse modeling and source detection. 

There is a need in the y-ray astronomy community for a method of separating source and 
background ‘diffuse’ emission in the Milky Way. In high-energy GeV and TeV data analysis, the 
challenge lies in finding methods that are suitable for use with the low statistics observed in this 
regime in order to provide reliable, stable results. 

Galactic emission is thought to be the sum of a number of distinct components [1]. However, 
for the purposes of this work we regard this as just two: the contribution from sources, and the back¬ 
ground flux. In this simplified picture, sources are defined as resolvable detections which fall above 
a given threshold, while the remaining flux (attributed to the contributions from unresolved sources 
and truly diffuse emission processes such as inverse Compton, decay and bremsstrahlung) forms 
the background - or ‘diffuse’ - emission. 

2. Method 

An algorithm based on convolution with an elon¬ 
gated background filter and significance clipping is in¬ 
troduced. Its operation is summarized in the following 
steps. 

• Background: Convolution of a counts map with 
a background kernel to act as an initial diffuse 
background estimation. 

• Significance: Computation of a significance im¬ 
age from the counts map and background image 
using the Li & Ma method [2]. 

• Mask image computation: Determination of 
an exclusion mask from the significance image. 

This removes regions above a set significance 
threshold which are replaced with a diffuse es¬ 
timation based on the flux around the excluded 
source within a background convolution kernel. 

The initial diffuse estimation is taken as the to- Figure 1: Iterative background method 
tal flux contribution, as determined from the provided demonstration, 
counts image. The significance and mask image com¬ 
putation steps repeat, with excluded regions dilated in each iteration for which they fall above the 
significance threshold. Once there is no change in the exclusion mask, the final diffuse estimation 
image is returned. This iteration process is illustrated by Figure 1, where contours indicate the 
shape of the exclusion mask. 

Six parameters in Table 2, p. 5 (the values indicated refer to the example application to an 
dataset discussed in sections 3 and 4) specify the method, which can be considered a generalization 
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of the ring background method [3] adapted to handle elongated diffuse estimation in the Galactic 
plane through the use of an lengthened background kernel. 

3. Dataset 

The datasets upon which the method is demonstrated comprise both of experimental observa¬ 
tions and simulated inputs. These are intended to provide a realistic and well-understood test-base. 
Both are described in this section. 

For all datasets, an analysis region of Galactic latitude —5° <b <5° and longitude —100° < 
/ < 100° is used to encompass the majority of the observed Galactic emission. 

3.1 Observational Data 

Fermi-LKE data from 5 years of observation was chosen with a well defined 10-500 GeV 
energy cut as a showcase for the method. This choice was due to a well-understood PSF, good 
coverage of the Galactic plane and roughly uniform exposure of the region of interest. This region 
and energy band contains 6.2 x 10"^ events, following a power-law of index T = —1.55. It should 
be noted that applications of this method in analyses need not be limited to Fermi-LP^E data: image 
analysis of all high-energy Galactic data is the intended use-case. 

3.2 Simulated Data 

The simulated dataset may be split twofold: inputs for which the source population was based 
on observation and inputs for which a modeled source population was used. In both cases, the 
sources were added to a known diffuse contribution to build a complete flux image. This was 
taken as the most recent Fermi diffuse background model gll_iem_v0 5_revl. fit, offering 
a realistic distribution of background flux throughout the full study region. 

Three different source populations were used. The first of these was a PSF-convolved point 
source image of the IFHL catalog above 10 GeV. The resulting source-only image had a flux 
of 5.06 X 10 ^ph cm ^s ^ within the study region. The remaining two source populations were 
entirely modeled, based on the study undertaken in the simulation work of the IFHL Catalog paper 
[4, p.59] which follows the essence of the study presented in [5]. This considered a reference galaxy 
of a similar source population and distribution to the Milky Way, and two further simulations with 
different densities and proportions of bright and faint sources. Here we considered the model 
of similar parameters to those measured in the Milky Way and a further case where simulated 
parameters were chosen to produce a model of similar flux distribution to the IFHL catalog above 
its detection threshold. 

In line with [5], it was assumed that the luminosity function of the y-ray source population 
between two limits Ly^min and Ly^max followed a simple power law. Further, it was taken that the 
spatial distribution of pulsars offered a representative sample of the full spatial distribution of a 
Galactic y-ray sources [5, p.2]. Thus a galactocentric {R^z) pulsar distribution model was employed 
[6, p.7] (a Gamma function) to distribute sources for the p{R) radial density distribution function. 
An exponential function was used for the p (z) distribution. Standard Monte-Carlo methods were 
used to sample p(Ly,/^,z) throughout the simulated galaxy. The resulting distribution was then 
normalized and scaled at Rq to the observed population density at the Sun. In line with the IFHL 
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Source Model 

Population Density 
at the Sun/kpc~^ 

Minimum Luminosity/ 
lO^^ph s“^ 

Maximum Luminosity/ 
10“ph s' 

A 

3 

10 

10 

B 

10 

4.0 

4.0 


Table 1: Parameters for 10 - 500 GeV Galaxy Population Simulations, 
study, the radial distribution was adjusted to peak at 4 kpc and the p(z) exponential scale was set 
to 0.5 kpc. 

A reference model (A) was defined with Pq = 3kpc“^, minimum y luminosity Ly^min = 
lO^^ph s“^ and maximum y-ray luminosity Ly^max — lO^^ph s“^ with luminosity power-law in¬ 
dex r = —1.5. The comparison model, B - of higher population density at the position of the Sun 
and lower luminosity bounds - was also considered. These are summarized in Table 1. 

Figure 2 shows the flux distribution of each of the 
input galaxy populations. In the case of the IFHL cata¬ 
log, the detection threshold of order 10“^^ph cm“^s“^ 
is evident. No such threshold is applied to populations 
A and B. Thus it can be anticipated that the separation 
algorithm returns better results in the case of the IFHL 
source population, as there will be fewer unresolvable 
faint sources. 

Before application to the source separation al¬ 
gorithm, the IFHL and simulated catalogs are added 
to the integrated Fermi Diffuse background model 
gll_iem_v0 5_revl. f it between 10 GeV and 500 GeV to produce the input source & diffuse 
test cases. The integral flux of this model in the Galactic plane region was 5.32 x 10“^ph cm“^s“^ 
The fluxes of the three source populations are included in Table 3. 

4. Results 

The diffuse estimation algorithm was applied to the described inputs with the parameter values 
presented in Table 2. The significance threshold was set according to preliminary studies so as to 
avoid significance-triggering (and hence false source detection) on Poisson up-fluctuations. The 
correlation radius and mask dilation radius were chosen to be of the order of the Fermi-hAI PSF to 
optimize for point source detection which would be observed on these scales (the LAT PSF has a 
68% containment radius of 0.150°and 95% containment radius of 0.867°in the 10-500 GeV energy 
band), while an elongated background kernel of suitable size was used to ensure all exclusion 
regions could be covered. A pixel size of 0.1°offered good resolution for source detection without 
leading to an excessive computational load during analysis. 

The results from simulated galaxy models are presented in Table 3. The performance of the 
algorithm can be assessed by comparing the true source flux fraction with the recovered source flux 
fraction. In the case of the IFHL catalog sources more source flux is separated by the algorithm 
than was input, indicating some cases of source confusion. This is because the diffuse model has 
small-scale structure and the process of smoothing some of the peaks in the background model 
leads them to be falsely detected as sources in the following iteration. This is particularly evident 



Figure 2: Simulated flux distributions. 
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Source Model 

Source Flux/ 
10“^ ph cm“^ 

Total Flux/ 

10“^ ph cm“^ 

True Source 

Flux Fraction/% 

Recovered Source 

Flux Fraction/% 

IFHL Catalog 

5.1 

58 

9 

16 

A 

8.5 

62 

14 

13 

B 

12.9 

66 

20 

10 


Table 3: Galactic Plane recovered diffuse fluxes 
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Figure 3: Background estimation for Fermi-hAI Data, 
around the galactic center. A reduction of this effect can be achieved by fine-tuning the algorithm 
parameters. For instance, raising the significance threshold would be expected to reduce incorrect 
source triggering due to diffuse background structure. 

It should be noted that there are instances in which this 
diffuse estimation algorithm will either fail or return incorrect 
results. Primarily, in cases for which the background box ker¬ 
nel is too small, there will be regions around sources where 
the exclusion mask dilates to a size larger than the box, lead¬ 
ing to inadequate data within the kernel from which a diffuse 
estimate can be computed. Additionally, in cases of very low 
counts (^(10) events in the background kernel) the results 
will be dominated by Poisson fluctuations. In these cases it may be found that the algorithm fails 
entirely if significance values cannot be calculated for the full analysis region, or a null result may 
be returned. 

For the simulated source populations, it is seen that source/diffuse separation deteriorates with 
as the proportion of fainter sources increases (model B has more faint sources). In these cases 
an effective threshold occurs where the significance of the sources against the background is not 
sufficient for detection. This threshold might be reduced through adjustment of the shape of the 
source kernel to better match the dimensions of PSF convolved faint sources, although in many 
cases the faintest of these sources will remain unresolvable by any method. 

Application to true Fermi-hAI data to develop a diffuse model (Figure 3) produced a diffuse 
estimation with a flux of 5.17 x 10“^ ph cm“^ s“^ in the Galactic region. Figure 4 presents spa¬ 
tial profiles of the resulting diffuse model. It can be seen that the diffuse estimation here largely 
follows the total flux distribution of the Milky Way with some sources excluded, notably Vela at 
-2.78°latitude and -96.4°longitude. Shortcomings are evident where diffuse estimation exceeds the 
total flux. This arises due to an inappropriately large size of the background convolution kernel in 
such regions, combined with source contamination in the diffuse estimation. 

5. Conclusions 

This work outlined a new method of modeling the Galactic diffuse emission and presented a 
diffuse emission model based on Fermi-LAT data between 10 and 500 GeV. This model is built 
only from the observed data, with the assumption that it is smooth on longitude and latitude scale 


Parameter 

Value 

Significance threshold 

4 G 

Correlation radius 

0.3° 

Mask dilation radius 

0.3° 

Height of background kernel 

0.5° 

Width of background kernel 

10° 

Pixel size 

0.1° 


Table 2: Algorithm Parameters 
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Figure 4: Fermi-LAI Profiles for total Galactic flux and estimated diffuse flux 
given by the background kernel. This is fundamentally different from the diffuse emission model 
usually used for Fermi-LAT data analysis [7], which is based on a much more complex physical 
model of cosmic rays and emission processes in the Milky Way [1]. 

Overestimation of the diffuse flux with this method was evident in some regions of resulting 
background estimations, and steps are suggested for fine-tuning. Combined with algorithm exten¬ 
sions to allow for the iterative use of multiple background kernels, improvements in the resulting 
diffuse estimation may be achieved. An extension of this could be to generalize the algorithm for 
multi-scale detection by computing significance images with multiple correlation radii such that 
source exclusion would occur for a wider range of morphologies, not just point like sources and 
those with a small angular size. The methods and algorithm introduced in this study, including 
the Galactic population simulation methods, are provided in the affiliated Astropy open-source 
package, Gammapy [8]. 
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